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ABSTRACT 


The basic structure of the solar system is set by the presence of low-mass 
terrestrial planets in its inner part and giant planets in its outer part. This is the 
result of the formation of a system of multiple embryos with approximately the 
mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer 
disk, within the lifetime of the gaseous component of the protoplanetary disk. 
What was the origin of this dichotomy in the mass distribution of embryos/cores? 
We show in this paper that the classic processes of runaway and oligarchic growth 
from a disk of planetesimals cannot explain this dichotomy, even if the original 
surface density of solids increased at the snowline. Instead, the accretion of drift¬ 
ing pebbles by embryos and cores can explain the dichotomy, provided that some 
assumptions hold true. We propose that the mass-flow of pebbles is two-times 
lower and the characteristic size of the pebbles is approximately ten times smaller 
within the snowline than beyond the snowline (respectively at heliocentric dis¬ 
tance r < Tice and r > rjce, where rice is the snowline heliocentric distance), dne 
to ice snblimation and the splitting of icy pebbles into a collection of chondrnle- 
size silicate grains. In this case, objects of original snb-lnnar mass would grow 
at drastically different rates in the two regions of the disk. Within the snow¬ 
line these bodies wonld reach approximately the mass of Mars while beyond the 
snowline they would grow to ~ 20 Earth masses. The results may change qnan- 
titatively with changes to the assnmed parameters, but the establishment of a 
clear dichotomy in the mass distribntion of protoplanets appears robust provided 
that there is enough turbulence in the disk to prevent the sedimentation of the 
silicate grains into a very thin layer. 



1. Introduction 


The Solar System has a characteristic structure, with low-mass rocky planets in its inner 
part, often called terrestrial planets, and giant planets (gas-dominated or ice-dominated) in 
the outer part. 

A census of protoplanetary disks in clusters with known ages shows that the dust 
emission (usually assumed to trace the abundance of gas) disappears in a few My (Haisch et 
al., 2001); this is also the timescale on which the emission lines diagnostic of gas accretion 
onto the central star fade away (Hartmann et ah, 1998). The fact that no primitive 
chondrite parent bodies seem to have accreted beyond 3-4 My (Kleine et al., 2005) suggests 
that the proto-Solar-System disk was not of exceptional longevity. 

Clearly, the giant planets had to form within the lifetime of the gas-disk because they 
accreted substantial amounts of hydrogen and helium (this is true also for Uranus and 
Neptune). Thus giant planets should have formed within a few My only. The commonly 
accepted scenario for giant-planet formation is the core-accretion model (Pollack et ah, 
1996). In short, a massive solid core accretes hrst and then it captures a massive atmosphere 
of H and He from the protoplanetary disk. The mass of all giant planet cores but Jupiter 
is around 10 Earth masses (Mq) (Guillot, 2005). Jupiter might have no core today 
(Nettelmann et al., 2008). However, there are several tens of Earth masses of “metals” 
(molecules heavier than H and He) in Jupiter (Guillot, 2005) and it is possible that part or 
even most of its primordial core has been eroded and dissolved into the atmosphere (Guillot 
et al., 2004; Wilson and Militzer, 2012). 

An estimate of the mass of the core needed for the accretion of a massive atmosphere 
is also provided by models. It is generally considered, since the work of Pollack et ah, 
that the core needs to exceed ~ IOMq. More precisely, the critical mass for the runaway 
accretion of the atmosphere depends on the rate at which the core accretes solids, on the 
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molecular weight of the atmosphere (Ikoma et ah, 2000; Hori and Ikoma, 2011) and the 
dust opacity in the envelope (Mizuno, 1980, Stevenson, 1982), which remains poorly known 
despite modern attempts to estimate the dust opacity self-consistently by modeling the 
aggregation of infalling grains (Movshovitz and Podolak, 2008; Ormel, 2014). The solid 
accretion rate could not be arbitrarily small, otherwise the core would not have formed in 
hrst place within the lifetime of the disk. This gives a constraint on the minimal mass of 
the core. Lambrechts et al. (2014; see their Fig. 7) showed that the core should have had a 
mass of at least 10 M 0 , if the ratio H 2 O/H 2 in the atmosphere was less than 0.6. Uranus 
and Neptune provide an indirect conhrmation of this estimate as a lower-bound for the core 
mass. In fact, they have a core of about 10-15 Mq (Guillot, 2005) and only a few Earth 
masses of H and He, which means that they either did not start runaway accretion of gas, 
or did so only at the very end of the lifetime of the disk. 

The situation for the terrestrial planets is completely different. There is a general 
consensus that the terrestrial planets formed from a system of planetary embryos and 
planetesimals (see Morbidelli et ah, 2012 for a review), although the details of how this 
happened can differ from one model to the other (Chambers and Wetherill, 1998; Chambers, 
2001; Agnor et ah, 1999; Raymond et ah, 2004, 2006; O’brien et ah, 2006; Hansen, 2009; 
Walsh et ah, 2011; Jacobson and Morbidelli, 2014). According to these models and to the 
interpretation of isotopic chronometers for terrestrial and lunar samples (Yin et ah, 2002; 
Jacobsen, 2005; Touboul et ah, 2007; Allegre et ah, 2008; Halliday, 2008; Taylor et ah, 
2009) the Earth took several tens of My to complete its formation, with a preferred timing 
for the Moon forming event around 100 My (Jacobson et ah, 2014). The minimum time 
in which the Earth acquired 63% of its mass is 11 My (Yin et ah, 2002; Jacobsen, 2005). 
Thus, most of the assemblage of the Earth clearly took place after the removal of the gas 
from the protoplanetary disk. 



5 


Mars, instead, formed very quickly, i.e. in a few My (Halliday and Kleine, 2006; 
Dauphas and Pourmand, 2011), basically on the same timescale of chondritic parent bodies. 
This suggests that Mars is a stranded embryo (Jacobson and Morbidelli, 2014). The fact 
that the Moon-forming projectile also had a mass of the order of a Mars-mass (Canup and 
Asphaug, 2001; Cuk and Stewart, 2012) and that Mercury, if it had originally the same 
iron content as the other terrestrial planets, was also approximately Mars-mass (Benz et 
ah, 1988) suggests that the mass of Mars was the typical mass of planetary embryos in the 
inner Solar System at the time the gas was removed from the protoplanetary disk. 

In summary, it appears compelling that, by the time gas was removed from the system, 
the process of formation of the solid component of planets had produced a great dichotomy 
in the mass distribution of protoplanets: in the inner system, the largest objects were 
approximately Mars-mass; in the outer solar system they were ~ IOMq. Thus, there was a 
contrast of two orders of magnitude between the masses of the solid planets formed in the 
inner and outer systems respectively. This happened despite the accretion timescale, which 
can be reasonably approximated by the orbital timescale, is 10 times faster at 1 AU than 5 
AU! 

How could this be possible? The generic (and hand-waving) explanation is that 
the cores of the giant planets formed beyond the ice line, so that the density of solids 
was comparatively larger than in the inner solar system, where only refractory material 
could be in solid form. However, this can not be the explanation. According to Lodders 
(2003), for the solar abundance the H20-ice/rock ratio is approximately one-to-one. That 
means that the amount of solid mass available for planet formation beyond the snowline 
increases by just a factor of 2. This is conhrmed by the ice/rock ratio inferred for comets, 
trans-Neptunian objects, and irregular satellites of giant planets (McDonnell et al. 1987; 
Stern, 1997; Johnson and Lunine, 2005). An enhancement of the solid mass by more than a 
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factor of 2 might have been produced by the so-called “cold-huger effect” (Morhll and Volk, 
1984; Ros and Johansen, 2013), but this would have happened only locally at the snowline 
and therefore could explain at most the formation of one giant-planet core, not several. 

The purpose of this paper is to investigate which process of planet formation is more 
likely to have led to the dichotomy discussed above. In section 2 we consider the classical 
process of formation of embryos/cores by runaway/oligarchic accretion of planetesimals 
(Greenberg et ah, 1978; Kokubo and Ida, 1998; Wetherill and Stewart, 1993; Weidenschilling 
et ah, 1997). We show that this process clearly cannot explain the dichotomy. Next, in 
section 3, we consider the process of pebble accretion. This is a new process for planet 
growth, introduced in Lambrechts and Johansen (2012; see also the precursor work by 
Ormel and Klahr, 2010; Johansen and Lacerda, 2010; Murray-Clay et ah, 2011; Bromley 
and Kenyon, 2011), which is rapidly gaining attention (Morbidelli and Nesvorny, 2012; 
Chambers, 2014; Lambrechts and Johansen, 2014; Lambrechts et ah, 2014; Guillot et ah, 
2014; Kretke and Levison, 2014a,b). We will show that, provided some assumptions hold 
true, the pebble accretion process can explain the two orders of magnitude mass-contrast 
between inner solar system objects and outer solar system objects. The conclusions and 
perspectives are discussed in section 4. 


2. Growth of embryos and cores by planetesimal accretion 

The growth of embryos and cores from a disk of planetesimals proceeds in two phases. 

The first phase is that of runaway growth (Greenberg et ah, 1978). Here most of the 
mass of the disk is in “small” planetesimals. The velocity dispersion of the planetesimals 
is set by the equilibrium between the self-excitation of their orbits, also called self-stirring, 
and gas drag. We neglect here collisional damping because it is important only for very 
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small objects, which we will call pebbles in the next section, and in the absence of gas 
drag (Goldreich et ah, 2004; Levison and Morbidelli, 2007). The velocity dispersion of the 
planetesimals is therefore comparable or smaller (because of the drag) to the escape velocity 
from the surface of the planetesimals carrying the bulk of the population mass (Greenberg 
et ah, 1978). In this situation, the accretion cross section a of an individual planetesimal is: 



( 1 ) 


where R is the planetesimal radius, Vesc is the escape velocity from the planetesimal surface, 
Vrei is the dispersion velocity in the planetesimal disk and the term in parenthesis is called 
the “gravitational focusing factor” (Greenberg et ah, 1978; Greenzweig and Lissauer, 1990, 
1992). Thus, the most massive planetesimals have a comparative advantage. If their 14sc 
is significantly larger than Vrei, their gravitational focusing factor can be approximated by 
Vesd^rei ^ where M is their mass. Because the accretion rate is proportional to 

a and R oc from ([^ the relative mass accretion rate of the massive bodies is: 


Eq. (|^ means that the most massive bodies grow the fastest and their mass ratio with 
the rest of the planetesimal population increases exponentially with time. Hence the name 
“runaway growth”. 

The second phase is that of oligarchic growth (Kokubo and Ida, 1998). It occurs 
when the runaway bodies become massive enough to stir-up the velocity dispersion of the 
planetesimals to a value comparable to the escape speed of the runaway bodies. In this case 
the gravitational focusing factor becomes of order unity and therefore eq. ([^ becomes: 



( 3 ) 


This equation means that the more massive a body is, the slower (in relative terms) it 
grows. The difference in mass between bodies becomes smaller and smaller as time passes. 





In principle, the planetesimals could reduce the mass difference with the former runaway 
bodies. But in general, when (|^ holds, the relative velocities are so large that collisions 
among the planetesimals become disruptive. Thus, only the massive bodies keep growing, 
forming a population of oligarchs that dominate the dynamical evolution in the disk. Hence 
the name “oligarchic growth”. Numerical simulations (Kokubo and Ida, 2000) show that 
the oligarchs accommodate their orbits (possibly by merging with each other) so that their 
separation in semi major axis is 5-10 mutual Hill radii. Oligarchic growth stops when an 
oligarch has accreted most of the mass available in the annulus around its orbit that extends 
up to half-way to the neighboring oligarch. 

In order to compare the masses of protoplanets growing in different parts of the 
protoplanetary disk, below we first compare the initial runaway growth rates, then the final 
masses at the end of the oligarchic growth phase. 


2.1. Radial dependence of the runaway growth rate 

Here we assume that, at the beginning of the growth phase, both the planetesimals 
and the runaway bodies have the same mass in all parts of the disk. This is in line with the 
classic assumptions that planetesimals are km-size everywhere in the disk and the runaway 
bodies are initially barely bigger than the characteristic planetesimal size (e.g. Kokubo and 
Ida, 1996; Greenberg et ah, 1984). Under this assumption, on which we will return at the 
end of this section, the radial dependence of the growth rate of the runaway bodies is set 
by the radial dependencies of the spatial density of the planetesimal population and of the 
orbital properties (i.e. velocity dispersion). 

First, let’s analyze how the spatial density of the planetesimal population pp depends 
on r. For this purpose we start from the radial profile of the gas density. 
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Whatever its initial distribution, the inner part of a protoplanetary disk of gas should 
rapidly evolve, under the action of viscosity, towards a distribution such that the mass-flow 
towards the star is independent of r (Lynden-Bell and Pringle, 1974). This is called an 
“accretion disk”. In this situation, and assuming the usual a-prescription for the viscosity 
(Shakura and Sunyaev, 1973), the exponent a of the power-law describing the density 
prohle of the disk oc l/r“ is: a = — 3/2 (Bitsch et ah, 2014a). Here Hg is the 

pressure scale height of the gas. Bitsch et ah (2014a, 2015) showed that the outer part of 
the disk is usually flared under the effect of stellar irradiation, so that =9/7 and 

a = 15/14. However, the inner part of the disk, dominated by viscous heating, is not flared. 
Thus ~ 1 (apart from some oscillations due to opacity transitions, which we neglect 

here). For an early disk, with a mass flow to the star of Mq/ij (where Mq is the mass 
of the Sun), the transition between these two regimes occurs at ~ 6 AU. Thus, for our 
purposes we consider that the disk is not flared, so that a ~ 1/2. The volume density of the 
gas Pg oc Tig/Hg F therefore proportional to 1/r^/^. 

This disk’s density prohle is much shallower than that in the usual Minimum Mass 
Solar Nebula (MMSN) model (Weidenschilling, 1977; Hayashi, 1981), in which E^ oc 1 /r^A 
and Pg oc The fundamental assumptions of the MMSN model (the planets formed 

in situ, exclusively from local material), however, are not considered valid any more (e.g. 
Crida, 2009) given what we now today know about planet migration; thus we think that the 
shallow disk prohle with a = 1/2 has a stronger theoretical motivation. It is also important 
to realize that if the disk’s density prohle is shallower, then the chances to form massive 
planets further out increases (e.g. Lissauer, 1987). Thus, using a shallower disk is more 
favorable for explaining the solar system dichotomy. 

At the beginning of the accretion phase, it is commonly assumed that the surface 
density of solids Ep is proportional to that of the gas: Ep = mE^, where m is the metallicity 
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factor. As we said in the introduction m should increase by a factor of ~ 2 at the snowline 
because of the presence of water ice (Lodders, 2003). The volume density of planetesimals 
on the disk’s midplane, Pp, is proportional to Sp/r sin(i), where i is the inclination dispersion 
of the planetesimals. It is well-known that for a disk undergoing self-stirring, i oc e/2 
(Stewart and Ida, 2000). Thus we have: 


Pp oc 




(4) 


We now analyze the velocity dispersion of the planetesimal disk, namely its eccentricity 
excitation. 


As we said above, the eccentricity excitation is set by the equilibrium between 
self-stirring of the planetesimal population and gas drag. For a disk of equal mass particles, 
the growth rate of the eccentricity dispersion is (Stewart and Ida, 2000; Morbidelli et al., 
2009 - supplementary material): 


de^ PpM 
df V^e 

where M is the individual mass of the planetesimals and Vk is the keplerian velocity. 


(5) 


For the gas drag, the evolution of the eccentricity excitation is (Wetherill and Stewart, 
1989; Morbidelli et al., 2009 - supplementary material): 


de2 epgV^R^ 
- (X - - — 

df VkM 


( 6 ) 


where R is the radius of the planetesimals and Vg is the orbital velocity of the gas which, 
for eccentric planetesimal orbits, can be safely approximated by Vk- 


Thus, by equating ([^ and (|^ and using the relationship (|^ and keeping only the 
terms which have a radial dependence, we hnd that the velocity dispersion is: 


Vrei = eVx oc 


K 


V, 


1/3 


oc r 


1/6 


K 


(7) 
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so that e oc 

We can now use these results in the accretion formula: 

M = PpaVrel , (8) 

with a given in Q. We hud: 

M oc —— oc . (9) 

eVrel 

This accretion rate is the same as that found in eq. 7 of Thommes et ah (2003), if one 
remembers that in that paper is the same as rpp here and one plugs in the radial 

dependences of Pp and e given above. 

The conclusion is that accretion is much faster in the inner disk than in the outer disk, 
even in a disk with the shallowest possible density radial prohle like the one we considered. 
With the radial dependence given in ([^, even taking into account that m is twice larger 
beyond the snowline, the accretion rate at 1 AU is 20 times faster than at 5 AU! Because 
in runaway growth the relative accretion rate M/M oc absorbing this factor of 20 

between the accretion rates at the two locations would require that the runaway body at 
5 AU is initially 8,000 times more massive than at 1 AU. We don’t see any justihcation of 
why it should be so. 


2.2. Radial dependence of the final mass of the oligarchs 

The completion of oligarchic growth is quite an idealized concept. Oligarchic growth is 
rather slow and therefore it takes time to bring the process to completion, easily exceeding 
the lifetime of the gas in the protoplanetary disk. The simulations of Levison et ah (2012), 
the most realistic to date as they combine embryos growth with planetesimal grinding 
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semi major axis (AU) 


Fig. 1.— The mass distribution of embryos predicted by (10), for parameters such that the embryo 


at 1 AU has a 0.1 Mm. The snowiine is set here at 3.5 AU. 


and account for a full size-distribution of objects, show that oligarchic growth approaches 
completion only in the very inner part of the disk. This result is in agreement with the 
strong radial dependence of the runaway growth rate found in the previous section. 

Nevertheless, it is instructive to compute which kind of radial mass distribution 
of objects oligarchic growth would predict, if it reached completion. By dehnition of 
“completion”, the mass of the body is of the order of that originally available in an annulus 
around the orbit with a width equal to n times the Hill radius Rh of the body itself, with 
n ~ 10, independent of heliocentric distance (Kokubo and Ida, 2000). Thus we can write 


Mali = 27rrni?HS„. 


( 10 ) 


With the shallow Sp oc ml^ dehned above and remembering that Rh = r{M^ 


eq. (10) gives Mou cx: mr^!^. 


This is a strong growing function of r. If one accounts for an increase of m by ~ 2 
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beyond the snowline, Moh predicts a ratio of almost 100 between the oligarchic masses at 
5 AU and 1 AU! Is this the sign that oligarchic growth reached completion thronghont the 
Solar System? 

We donbt that this is the answer for two reasons. First, look at Fig. which shows 
the mass distribntion of the oligarchs predicted by the formnla above. Here, the snowline 
effect has been assnmed at 3.5 AU. The formnla does not predict Mars-mass embryos 
through the inner solar system and massive cores beyond the snowline. It predicts a mass 
distribution smoothly growing with distance up to the snowline, with Earth mass bodies in 
the asteroid belt. We don’t see how this distribution of bodies could be reconciled with the 
current structure of the Solar System, given current models of terrestrial planet formation, 
planet migration etc. Second, Levison et ah (2010) showed that when bodies achieve a few 
Earth masses in the outer part of the disk their growth slows down enormously, because 
they start to scatter away the planetesimals from their feeding zone, rather than accreting 
them. Thus, all masses predicted in Fig. [^beyond 3 AU are most likely overestimated and 
~ IOM 0 cores are unable to form. 

In summary, we conclude that the classical scheme of runaway/oligarchic growth of 
protoplanets from a disk of planetesimals can not explain the great dichotomy of the Solar 
System and hence its current structure. In the next section we turn to a new paradigm for 
planetary growth, that of pebble accretion (Lambrechts and Johansen, 2012). 


3. Growth of embryos and cores by pebble accretion 

In a series of recent papers (Lambrechts & Johansen, 2012, 2014; Ormel and Klahr, 
2010; Johansen and Lacerda, 2010; Murray-CIay et ah, 2011; Morbidelli and Nesvorny, 
2012; Kretke and Levison, 2014b), a new model has been proposed in which massive 
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planetesimals accrete pebbles (particles cm to dm in size) very efficiently, thanks to the 
combination of gravitational deflection and gas drag. If most of the solid mass of the disk 
is in pebbles, the largest planetesimals can rapidly grow to several Earth masses. 


The process of pebble accretion can be illustrated and quantihed as follows. First, 
it is convenient to dehne the Bondi radius Rb as the distance at which a planetesimal of 
mass M, assumed to be on a circular orbit, exerts a deflection of one radian on a particle 
approaching with a speed Av equal to the relative speed of the ga^ . Because the gas is 
pressure supported, its orbital velocity is sub-Keplerian, and therefore it approaches the 
planetesimal with a speed An = tjVk with: 

- _1 dlogP 

^ 2 \ r ) d log r ’ 

where P is the gas pressure. Denoting by G the gravitational constant, one has: 



Rb — 


GM 
(An)2 • 


( 12 ) 


Due to three-body effects, however, the deflection cannot be effective beyond the Hill radius 
Rh = r[M/(3M0)]^/^ of the planetesimal. Thus, we dehne : 


Rgp — min(i?B, Rh) , 


as the distance from the planetesimal within which its gravitational pull is effective. 
Because the Bondi radius is proportional to M while the Hill radius is proportional to 
for small bodies Rgp = Rb, while for more massive bodies Rgp = Rh- For a typical value 
of rj of 0.0027-0.0032, one has Rb = Rh for M ~ 2-3xlO“^M0. Here we will consider for 

^Notice that this dehnition of the Bondi radius, given in Lambrechts and Johansen (2012) 
is different from the classic dehnition of Bondi radius given in studies of bound planetary 
atmospheres, namely the distance from the planet at which the escape speed is equal to the 
sound speed. 
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simplicity only bodies more massive than this threshold, so we detail the process of pebble 
accretion only in the case Rgp = Rr- 

Lambrechts and Johansen (2012) find analytically and numerically that for pebbles 
with Stokes number r (the product between the friction time -the timescale on which a 
particle looses by friction its velocity relative to the gas- and the orbital frequency) the 
effective radius for accretion onto the planetesimal is: 

reff = (r/0.1)^/^i?GP , (13) 


which is valid for r < 0.1. Notice also that the formula above has been derived for accreting 
bodies on circular orbits, for which the velocity difference with the pebbles is due to the 
Keplerian shear. If the body is eccentric, its velocity relative to the circular orbit has to be 
considered in the calculation of the relative velocity of the pebbles and, when the relative 


velocity increases, the value of re// drops. In Sect. 3.5.2 we will confirm that the eccentricity 
excitation of accreting embryos can be neglected. 


To compute the growth rate on a body, one needs to distinguish between 2D and 3D 
accretion, depending on the pebbles scale height Hp^ relative to rg//. Below we review the 
basic formulae for 2D and 3D accretion, then decide what is the threshold for Hp^ to pass 
from one regime to the other one. 


2D accretion 

Let’s consider that Hph is small relative to rg//. The accretion rate of the body is then: 

^2D ‘^'^eff^rel^pb 

where is the surface density of the disk of pebbles and Vrei is the relative velocity 
between the pebble and the body, namely 


Vrel = max{Av,V shear) 


( 14 ) 



16 


and Vshear = ^ = a/GMqT^- Because pebbles undergo inward radial drift due 

to gas drag, there is a continuous mass flux Mp of pebbles through any orbital radius. 
Thus, the steady-state surface density of pebbles is related to this mass radial flux via the 
relationship: 

Spfe = Mp/{2nrVr) 

where is the radial drift velocity of pebbles, which can be approximated for small r as 
Vr = 2tAv. 


Thus we have: 


M 2 D — 


2r effVrelMp 

AttttAv 


(15) 


3D accretion 

In this case we have: 

Msd = nrlffVrelPpb (16) 

where Pph is the volume density of pebbles on the midplane, which is related to the surface 
density by the relationship: 

Ppb Bpb/( v^27rh/p^) . 


Regrouping terms we hnd: 


Msd — M2D 


f \ 


Transition from the 2D to the 3D regime 

From the last formula, it is obvious that the transition occurs when: 

^reff 

2\/^Hpb 


( 17 ) 


< 1 . 
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This sets a condition on Hpj,. The latter is related to the disk scale height Hg by the 
relationship: 

Hpb = Hg^/(^ , (18) 

(Youdin and Lithwick, 2007), where a is the turbulence parameter of the disk of gas in the 
Shakura and Sunyaev (1973) prescription. 

Thus, 2D accretion will occur only for large planets (i.e. large r^ff) in disks with thin 
particle layers (i.e. small Hpb). This is the case if turbulence is weak (i.e. small a) and 
particles large (i.e. large r). For small planets and thick particle disks, accretion is 3D. 


Dependence of the accretion rate on r 


We now develop the dependence on the particles’ Stokes number r of the accretion 
rates reported in Eq. (15) and For this purpose, we consider Erst the most likely case 


where Vrei = vshear (see eq. 14). For the case of 2D accretion, because r^ff oc (eq. 13) 


and Vshear OC Tg// (eq. 14), formula (15) gives: 

M 2 D oc . 


For the case of 3D accretion, because Hpb oc r (eq 18), formula (17) gives: 

oc . 


(19) 


( 20 ) 


If instead Vrei = An (this happens only when rg// is very small), it is easy to see that the 
expressions for M 2 D and M^p become proportional to ^nd respectively. 


Thus, it is important to realize that in 3D accretion, the smaller is the value of r (i.e. 
smaller particles or higher gas density) the less efficient is the accretion. In 2D accretion, it 
is the opposite. 


With these formulae in hand, we can now design a scenario which may explain the 










great Solar System dichotomy and test it with simulations. 


3.1. Model scenario and simnlations 

Planetary embryos and cores grow from the pebbles that flow through their orbits. 
The basic assumption of the model we propose is that the pebble flux is quite different 
within and beyond the snowline (respectively for r < rice and r > r^ce where r^ce is the 
location of the snowline). Beyond the snowline there is a flow of icy pebbles. Lambrechts 
and Johansen (2014) have developed a pebble coagulation and drift model based on the 
work of Birnstiel et ah (2012). They found that, at t = 10^-10® y of the disk’s evolution, 
icy pebbles at ~ 4 AU are a few cm in radius with r ~ 10“^'®. They carry a mass flux of 
Mp ~ Me/ y. Here we will assume that the snowline is at 3.5 AU, consistent with the 
Grand Tack model of evolution of the giant planets and formation of the terrestrial planets 
(Walsh et ah, 2011; Jacobson and Morbidelli, 2014), but the exact location of this line does 
not have practical importance in what follows. 

When the pebbles come to the snowline, it is unlikely that they stop drifting towards 
the Sun. The drift of particles can stop only in presence of a pressure trap, where the radial 
pressure gradient of the gas is positive, so that the gas is locally super-Keplerian. It has 
been proposed (Kretke and Lin, 2007) that the snowline can create a pressure bump because 
the condensation of icy grain could strongly reduce the viscosity of the gas, thus increasing 
the gas surface density at the snowline location. However, improved 3D simulations of the 
disk structure in presence of viscosity transitions (Bitsch et ah, 2014b) show that this is 
true only in presence of extreme viscosity transitions, with < —28 where u is the 

viscosity. It is very unlikely that the condensation of icy grains can have such a dramatic 
effect on the viscosity. 
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Thus, we assume that the icy pebbles drift past the snowline. When they penetrate 
into the warm part of the disk, the ice sublimates away. This decreases the muss flux by 
a factor of ~ 2 (Lodders, 2003). Moreover, we assume that the icy pebbles are a mixture 
of ice and numerous silicate grains, hold together by the ice; thus, by sublimation, the icy 
pebbles release silicate grains which are much smaller than the original icy pebbles; these 
grains have a signihcantly smaller Stokes number and are much more coupled to the gas. 

It is unclear a priori what the typical size of the silicate grains could be. However, 
observations can give us some hints. Ordinary chondrites are made mostly of chondrules, 
which are mm-sized particles (Friedrich et ah, 2014). We do not know whether the 
sublimation of the icy pebbles released chondrules, or just the precursor grains of 
chondrules, but in terms of size the two are probably very similar. Moreover, coagulation 
experiments show the effectiveness of a mm-bouncing barrier for silicate particles (Guttler, 
2009), so it is likely that silicate particles could grow to mm-size relatively easily but not 
beyond (unlike icy particles which can easily grow to larger sizes; Wada et ah, 2009; Ros 
and Johansen, 2013). For all these reasons, we assume that the silicate particles released 
by the sublimation of the icy pebbles are millimeter in size. These particles would have a 
Stokes number r ~ 10“^'® at 3.5 AU. For simplicity we assume a unique size (i.e. value of 
r) for the icy pebbles and the silicate grains, rather than a distribution peaked around the 
considered values. 

We consider the mass growths of two embryos on either side of the snowline, at 3.5 Aig 
Both bodies have initially the same mass, which we set to half lunar masses (O.OO5M0),or 
about to 2.5 Pluto masses. This choice of the initial mass is rather arbitrary, but motivated 

^ Obviously it makes no sense that both bodies are at 3.5 AU. We do this assumption 
to suppress all dependencies on heliocentric distance and highlight only the consequences of 
the snowline effects on pebble size and mass flux. 
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by a few considerations: (a) we wish the masses of the two bodies to be the same, so that 
the observed differences in the growth rate are dne solely to the snowline effects and not to 
a difference in the initial masses; (b) we wish a mass signihcantly smaller than the target 
masses for the embryo and the core of the giant planet (0.1 and ~ IOM 0 , respectively); (c) 
but at the same time we wish the bodies to be massive enough to be from the beginning in 
the Hill regime, to avoid the complication of the transition from the Bondi regime (for the 
typical value of rj = 0.0027 that we assume here, the transition between the two regimes 


occurs at M = O.OO 37 M 0 ). An additional motivation will be given in sect. 3.4 Modeling 
the formation of these bodies is beyond the scope of this paper. However, a recent work 
(Johansen et ah, 2015) showed that they can be produced by a combination of streaming 
instability and pebble accretion. Our paper focuses on how subsequent pebble accretion, 
with different mass fluxes and pebble sizes, can lead to the great Solar System dichotomy 
between giant planet cores and terrestrial planet embryos. 


Because the dependence of the accretion rate on r is different in the 2D and 3D 
accretion regimes, we consider separately two cases: hrst that of a turbulent disk, for 
which the particle layer is thick and most of the accretion occurs in 3D and then that 
of a low-turbulent disk, for which the particle layer is thin and accretion is 2D from the 
beginning. 


3.2. Turbulent disk 


Embryos with M = O.OO 5 M 0 at 3.5 AU have a Hill radius Rh = 0.006 AU. We assume 
an aspect ratio for the disk of gas of 0.05 and a turbulence parameter a = 10“^. Thus, the 


scale-height of the disk of icy pebbles is 0.03 AU (see eq. 18) and that of the silicate grains 


is ~ 3 times larger. So, both embryos start growing in the 3D regime. 
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Fig. 2.— Mass growth due to a flux of pebbles of two embryos located on opposite sides of the 
snowline, here assumed to be at 3.5 AU. The parameter rj is assumed to be 0.0027. The red curve 
is for an embryo beyond the snowline, with a flux of 2 x 10~^ME/y in icy pebbles with r = 10“^'®. 
The green curve is for the inner embryo, for which the mass flux is reduced to 1/2(1 — F) relative 
to the outer embryo, and the pebbles’ Stokes number is r = 10“^'^. Here F is the fraction of the 
pebble flux which is intercepted by the outer embryo. When the outer embryo reaches 2 OM 0 the 
flux of pebbles towards the inner embryo stops (Lambrechts et ah, 2014). Thus, the inner embryo 
stops growing, and the simulation is stopped at this point. According to Lambrechts et ah, 2014), 
the outer embryo then starts accreting gas efficiently, to become eventually a giant planet. 

The accretion history of the two embryos is shown in Fig. The accretion of the inner 
embryo is penalized relative to that of the onter body for three reasons. First, the inner 
body sees ronghly mm-sized particles, which correspond to a valne of r 10 times smaller 
than of the pebbles beyond the snowline. Remember that in 3D accretion oc 
with f3 = 1/6 or 1/2 depending on the the scaling of Vrei] so a smaller r implies a smaller 
accretion rate. Second, the total mass flnx of particles is assnmed to be two times smaller 
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due to ice sublimation. Third, the flux of pebbles inside the orbit of the outer embryo is 
(1 — F) times that going towards the orbit of the outer embryo, where F is the Altering 
factor (Morbidelli and Nesvorny, 2012), namely the fraction of the pebble flux intercepted 
by the outer embryo, which is computed numerically from the embryo’s accretion rate. 

Thus, it is no surprise to see in Fig. [^that the outer embryo grows much faster than 
the inner one. The ratio between the initial growth rates is ~ 3. During their accretion 
histories, the outer embryo passes into the 2D accretion regime at 1.5 My. The inner one 
remains in the 3D regime all the time. The simulation is stopped when the outer embryo 
reaches 2OM0. This is because embryos this massive can open a partial gap in the gas 
distribution, enough to reverse the pressure gradient at its outer border. This reversal in 
the pressure gradient makes the disk locally super-Keplerian and the inward radial drift 
of pebbles stops there (Lambrechts et ah, 2014). Consequently, the inward embryo has 
to stop accreting, because the pebble flux is interrupted. Here, the outer embryo reaches 
2 OM 0 at t = 1.75 My, but this time can be changed by changing the pebble mass flux (here 
assumed to be constant over time for simplicity and equal to 2 x Me/ y)- The time 
scales linearly with the inverse of the assumed mass flux. Thus, what is important is not 
the time at which the outer embryo reaches 2 OM 0 but the mass that the inner embryo has 
when this happens. As one can see, the final mass of the inner embryo (green) is of the 
order of a Mars mass. Thus, we have reproduced the great dichotomy of the Solar System. 

We have tested different simulation set-ups in order to understand the role played by 
each of our assumptions. Without reducing the value of r by an order of magnitude for the 
pebbles inside of the snowline, the inner embryo would have grown 3.5 times more massive. 
Without reducing the overall pebble mass flux by a factor of 2 at the snowline, the inner 
embryo would have grown five times more massive (the final mass does not scale linearly 
with the mass flux because a more massive object accretes more efficiently). Without both 
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these reductions (the only difference between the inner and outer embryos now being just 
the filtering factor 1 — F), the inner embryo would have grown 100 times more massive 
(i.e. would have reached 10 M^). Restoring the nominal set-up concerning the pebble flux, 
but reducing the initial masses of the embryos by one order of magnitud^ we find that 
the outer embryo would have reached 20 in 2.35 My, but the inner embryo would have 
grown only to 1.5 x 1O“^M0. If the initial mass of the embryos is this small, bringing the 
inner embryo to Mars-mass requires no reduction in r inside of the snowline. 


3.3. Quiescent disk 


Let’s now consider a disk with a very low tnrbulence, with a = 10 In this case, the 


height of the icy particle layer at 3.5 AU is 0.001 AU (see eq. 18) and that of the silicate 
grains is ~ 3 times larger. In this case accretion is 2D for both the embryos from the 
beginning of the simnlation. In this case, the inner embryo has a comparative advantage 
relative to the onter one, because M 2 D oc with {3 = —1/3 or —2/3 depending on the 
scaling of Vrei] thus smaller particles result is faster accretion. There is nevertheless a factor 
of 2 redaction in the total mass flow of solids, dne to ice snblimation. The filtering factor 
(1 — F) is negligible at the beginning of the simulation. 


Fig-i shows the result in this case. The inner embryo starts to accrete slightly faster 
than the onter one. It is only later, when the outer embryo exceeds the mass of Mars, that 
the factor (1 — F) starts to become non-negligible; this rednces the growth rate of the inner 
embryo and the mass ratio between the onter and the inner embryos starts to increase. 

At the end of the simulation, when the onter embryo reaches 2 OM 0 , both embryos are of 


^ Notice that these bodies initially would accrete in the Bondi regime, i.e. Rgp = Rb- 


Here, we still apply formnla (13) even if accurate only for the Hill regime. 
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Fig. 3.— The same as Fig. but for a much thinner particle disk. Here the height of the layer 
of icy particles at 3.5 AU is 0.001 AU, and that of silicate particles is ~ 3 times as much. This 
corresponds to a very low turbulence in the disk, with a = 10“®. The result fails to explain the 
great dichotomy of the solar system. 

comparable mass. Clearly, this case does not reproduce the great dichotomy of the Solar 
System. 

Therefore, our explanation is valid only if the proto-solar disk was not so quiet. A value 
of a = 10“® is very small, even for a so-called dead zone where the Magneto-Rotational 
Instability is not active (Stone et ah, 2000). In fact, some turbulence can be sustained in 
the dead zone due to a number of effects, e.g. diffusion of the magnetic accretion stress 
(Turner et ah, 2007), baroclynic instability (Klahr and Bodenheimer, 2003), vertical shear 
instability (Nelson et ah, 2013). In particular, Stoll and Kley (2014) estimated that the 
vertical shear instability can sustain turbulence with a strenght equivalent to a parameter 
a of about few times 10“^. 
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Setting a = 10“^, we find that the inner embryo grows only to twice the mass of Mars, 
with an accretion history starting (and remaining) in the 3D regime. For a smaller than 
this threshold, a significant part of the accretion history of the inner embryo becomes 2D 
and a clear difference in mass relative to the embryo beyond the snowline disappears, as in 
Fi.g[^ Thus, in the rest of this paper we will assume a = 10 with the understanding 
that the results are quite similar for a > 10“^. 


3.4. Multiple embryos in the inner Solar System 

The successful simulations of terrestrial planet accretion (e.g. O’brien et ah, 2006; 
Walsh et ah, 2011; Jacobson and Morbidelli, 2014) require that all embryos in the inner 
solar system had a mass of the order of the mass of Mars, not just the one that was the 
closest to the snowline. Thus, in this section we examine the radial dependence of the mass 
growth of a chain of embryos located from 0.7 AU up to the snowline location at 3.5 AU. 

Here, there are again a number of assumptions to make. We assume that initially 
there are 40 embryo “seeds” with a mass of 0.005 M^, equally spaced between the two 
distance limits reported above. We chose the number 40 because this is approximately the 
number of embryos used in simulations of terrestrial planet formation (e.g. Jacobson and 
Morbidelli, 2014). We chose the mass of 0.005 because it is generally accepted that the 
inner solar system contained about 1,000 Ceres size bodies (Wetherill, 1992) with a mass 
~ 2 X 1O“^M0 so that, assuming a cumulative size distribution N{> D) (x 1/D^ like that 
of the largest main belt asteroids, we expect 40 bodies with masses > O.OO5M0. Obviously 
we do not expect that all embryo seeds had exactly the same mass, but in this experiment 
it is important that all masses are equal at the beginning of the simulation so that the 
observed differences in growth rates are due solely to pebble accretion effects and not to 
initial mass differences. Moreover, as the silicate grains drift towards the Sun we assume 
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that they do not evolve in size by coagulation or fragmentation. This is because there is no 
clear evidence for the dependence of the chondrule’s size with the heliocentric distance of 
the chondritic parent bodies (Friedrich et al., 2014). However, because the density of gas 
changes, the Stokes number decreases as the particles decrease their heliocentric distance. 
Given the shallow radial density prohle of the gas dehned in section 2 we assume that the 
Stokes number of a particle is r oc ^/r. 

In a hrst step, in order to understand the role of each parameter, we assume that r] in 
( pT| is independent of radius (r/ = 0.0027 as in the previous simulations). From formula 
(0. one can see that the parameter r] would indeed be constant if a disk had Hg oc r and 
P were a power-law function of r. The results is shown in Fig. which reports in red 
the growth of the giant planet core beyond the snowline (same as in the previous hgures), 
in green the growth of the outermost embryo at 3.5 AU (same as in Fig. [^, in blue the 
embryo in the middle of the chain (at 2.1 AU) and in magenta the innermost embryo (at 
0.7 AU). As one can see, the growth rate decreases slightly approaching the Sun. This is 
due to the decrease in r, but also because each embryo receives a flux of pebbles reduced 
by the other embryos located upstream relative to its position. In other words, counting 
the embryos from the innermost one, the pebble flux across the orbit of embryo is 
= n?2+i(i - Fj), where Fj is the hltering factor of embryo 

However, in a realistic disk, rj is not fully independent of r, because the aspect ratio 
of the disk of gas Hg/r has oscillations as a consequence of the opacity changing with 
temperature (Bitsch et ah, 2014a, 2015). The same is true for the radial dependence of r. 
Fig. I shows theses dependence for the disk studied in Bitsch et al. (2015) corresponding to 
an accretion rate on the star of 3.5 x 10“®MQ/y. We chose this disk among those studied 
in that work because it has the snowline at 3.5 AU, as we assumed above. The hgure also 
reports various analytic hts, valid in specihc ranges of heliocentric distances, that we will 
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Fig. 4.— Top: the same as Fig. but now showing the evolution of 3 embryos in a chain of 40. 
The red curve is for the giant planet core and the green cnrve is for the embryo #40 (the outermost 
one at 3.5 AU), both already shown in Fig. The blue curve is for the embryo #20 at 2.1 AU, 
and the magenta curve is for the embryo #1 at 0.7 AU. Bottom: the final mass distribution of the 
embryos as a function of heliocentric distance. 




- 28 - 


ro 

I 

O 



heliocentric distance (AU) 



Fig. 5.— Top: the dots show the radial dependence of rj according to Bitsch et ah (2015) 
for a disk with M = 3.5 x Mq/j. The colored solid lines represent analytic fits valid over 
specific intervals of heliocentric distance. Fit 1: r] = 0.002487; Fit 2: rj = 0.002487y^r/1.7; Fit 3: 
rj = 0.00328/(r/3.25). Bottom: the same but for the radial dependence of r. Here we assume a 
jump in pebble size by a factor of 10 at 3.5 AU, as illustrated by the black curve connecting the 
dots. Here the fits for r < 3.5 AU are: Fit 1: r = 0.00133-y/r; Fit 2: r = 0.00093 r. 
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Fig. 6.— The same as Fig. but now assuming r/(r) and r(r) as given by the fits shown in Fig. 


use below. 


For rj we assume Fit 1 for r < 1.7 AU, Fit 2 for 1.7 < r < 3.2 AU and Fit 3 for 
3.2 < r < 4 AU. For r we assume Fit 1 for r < 2 AU and Fit 2 for 2 < r < 3.5 AU. The 
value of r for the core placed just beyond the snowline remains equal to 10“^'®. With this 
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set-up, the results are shown in Fig. The results are quite similar to those shown in 
Fig-i The mass distribution of the embryos has a clear kink at 3.15 AU, where the r] 
function has a local maximum. In fact, between 3.15 and 3.5 AU, rj sharply decays with 
increasing heliocentric distance r. Pebble accretion is strongly favored if rj is small. Thus, 
the mass distribution of embryos increases sharply with r. Instead, for r < 3.15 AU, r] 
increases with increasing r. In principle, this would give a mass distribution of embryos 
decreasing with increasing r. But remember that, as the outer embryos start to grow, their 
hltering factors Fj grow as well and the inner embryos receive a progressively reduced flux 
of pebbles. Thus, the hnal mass distribution is still a moderately growing function of r, but 
not as steep as beyond 3.15 AU. 


3.5. Caveats and Issues 

We discuss here some aspects of the problem which have been neglected so far: disk 
evolution, the effect of self-excitation of the orbits of the embryos, the effect of radial 
migration and the mass distribution of the giant planet cores. 


3.5.1. Disk evolution 

Over time the disk evolves and loses mass. The stellar accretion rate decreases 
(Hartmann et al., 1998). Consequently, viscous heating drops and the transition between 
the inner region of the disk with roughly constant aspect ratio and the outer flared region 
dominated by stellar irradiation moves towards the star (Bitsch et ah, 2015). This effect 
changes the dependence of r; as a function of radius. The snowline location also moves 
towards the star. Moreover, with decreasing gas density, r increases for a given pebble 
size but on the other hand only smaller pebbles may survive in the disk (Lambrechts and 
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Johansen, 2014). 

In this paper we have neglected all of this. The goal here was to give a proof of principle 
that pebble accretion can explain the dichotomy in the mass distribntion of embryos and 
cores within and beyond the snowline and for this reason we preferred to remain simple, 
working in the framework of a snapshot of the disk’s structnre, taken more or less in the 
middle of its evolntion. A more qnantitative modeling of planet accretion acconnting for 
disk’s evolntion remains to be done and is left for a fntnre work. 


3.5.2. Embryos’ self-excitation 


Kretke and Levison (2014b) pointed ont that if the embryos acquire orbits with 
signihcant eccentricities or inclinations they stop accreting pebbles. In fact, the relative 
velocities between the embryos and the pebbles become large, which drastically decreases 


the Bondi radius (12). They proposed that this process is the key to explain why only a 
handful of giant planets cores formed. In their simulations a system of multiple proto-cores, 
while growing, starts to self-excite their orbital distribution. Because of dynamical friction, 
the most massive proto-cores remain on the most circular and co-planar orbits, while 
the less massive ones acquire larger eccentricities and inclinations. This shuts down the 
accretion of the latter, allowing only a few cores to grow to large masses. 


Given the population of multiple embryos that we consider here in the inner solar 
system, it is necessary to verify whether their mutual self-excitation of the orbits could 
change drastically their accretion history. The simulations presented above are not 
dynamical simulations. They just implement the analytic formulae described in Sect. 
Thus, we address the issue of orbital self-excitation with a separate N-body simulation. 


We assume a system of 40 Mars-mass embryos, from 0.7 to 3.5 AU, thus separated by 
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6.75 mutual Hill radii (the mutual Hill radius being defined as (a + a')/2[(M + M')/(3M0)]^/^ 
where a, M, and a', M' are the semi major axes and masses of two adjacent embryos). 
Their initial orbits are circular, inclined by 10“^ radians and with random orientations (to 
avoid that the problem is planar). We don’t include Jupiter in this simulation because at 
the stage we are interested in Jupiter was just a core of less than 20 Earth masses. 


We simulate the evolution of this system of embryos for 2 My with a version of the 
Symba integrator (Duncan et ah, 1998) modified to account for the tidal damping of 
eccentricities and inclinations exerted by the embryos’ gravitational interaction with the 
disk of gas. For massive bodies like the embryos, this is the most important damping effect, 
whereas gas drag is negligible (the opposite is true for planetesimals). The tidal damping 
formulae are those issued in Cresswell et ah (2007) and Cresswell and Nelson (2008) from 
hydrodynamical simulations. For this simulation, we assume that disk of gas has a density 
of 2400g/cm^ at 1 AU, corresponding to the MMSN, but with a radial profile of I/a/t 


expected for an accretion disk (see sect. 2.1). We neglect the disk torque inducing radial 


migration (migration being discussed in sect. 3.5.3), although there is some radial migration 
associated with the eccentricity damping (for instance the innermost embryo is at the end 
of the simulation at 0.69 AU and the outermost one at 3.36 AU). 


The system of embryos remains on stable orbits. No merger events have been recorded. 
At the end of the simulation the embryo with the largest eccentricity has e = 0.0018. The 
median value of the eccentricities of the embryos is 0.0006. All inclinations are below 10“^ 
degrees. Thus, the orbital excitation in terms of deviation from the Keplerian velocity 
(eUx) is smaller (in most cases much smaller) than the differential velocity of the gas 
relative to the embryos {t]Vk), with rj given in Fig. [^. Keep in mind that this simulation 
likely overestimates the real orbital excitation of the embryos. In fact, during most of their 
growth, the embryos are less massive than one Mars-mass and, even at the end of the 
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simulation illustrated in Fig. not all embryos have reached this value. Thus we conclude 
that self-excitation of the embryo’s orbits is not an important issue affecting their growth. 
Self-excitation would have become important if the embryos had grown well beyond the 
mass of Mars, as it is the case beyond the snowline. Thus our results do not contradict 
those of Kretke and Levison (2014b). It is because growth by pebble accretion is strongly 
reduced within the snowline that the inner solar system ended up with numerous small-mass 
embryos, while the outer solar system, where pebble-accretion is more vigorous, produced 
large cores but only in small number. 


3.5.3. Migration 

In all the calculations presented in this paper, as well as in the N-body simulation, we 
have neglected planet migration. Planet migration is due to the tidal interaction with the 
disk of gas. In hrst approximation, the migration speed is proportional to the planet’s mass 
(Tanaka and Ward, 2002). For the gas disk considered to build Fig. a Mars-mass object 
would migrate from 1.8AU down to 1.3AU in ~ 3 Myr. A migration range of 0.5 AU looks 
large compared to the typical embryo-embryo distance of ~ 0.1 AU. However, remember 
that: (i) all embryos would migrate at a comparable rate, so their mutual distances would 
change much more slowly and (ii) embryos are much smaller than Mars-mass for most of the 
lifetime of the disk, thus they would migrate much less then reported above. In conclusion, 
we think we can safely neglect their migration during the growth process at the stage of a 
proof of concept paper as this one. 

The situation for the giant planet core is different. According to the classical Type-I 
migration model (Goldreich and Tremaine, 1980; Tanaka and Ward, 2002) the core should 
migrate into the inner solar system very quickly. But a new result by Benitez-Llambay et 
al. (2015) shows that a fast-accreting body receives a positive torque from the disk, due to 
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the fact that it is very hot and it heats the surrounding gas in an asymmetric way (given 
that the disk is sub-Keplerian). Thus it is possible that the core refrains from migrating 
inwards. Once the core becomes more massive than a few Earth masses, this “accretional 
torque” should lose relative importance. But at this stage, the thermal structure of the 
disk becomes relevant, as it can excite an entropy-driven corotation torque on the core 
(Paardekooper and Mellema, 2006; Paardekooper et ah, 2010, 2011; Masset and Casoli, 
2009, 2010). Bitsch et ah (2014a, 2015) showed that the core should position itself a few 
Alls beyond the snowline (approximately where rj reaches the minimum value; see Fig. |^), 
at a location where all the different torques that it suffers from the disk cancel each other 
out. 

As the disk evolves and cools, the mass range of cores for which the torques cancel 
out shrinks (Bitsch et ah, 2015) so that, eventually, a core is released to Type-I migration. 
Thus, the assumption that the core grows to 2 OM 0 without migrating is, by no doubt, 
simplistic. We have done it because migration is complicated (and not fully understood, 
hence the appearance of new effects as the “accretional torque”) and we wanted to keep 
the model simple, as a proof of concept. The issue of core migration will be the object of a 
future work. Nevertheless, we can not resist to make a few speculations about the effect of 
migration at the very end of this paper. 


3.5.4- Cores mass distribution 

In the present paper we have considered only one core beyond the snowline. But at 
least 4 cores formed in the outer solar system. For the goal of investigating the mass 
contrast between the embryos within the snowline and the hrst core outside of the snowline 
it was enough to £x the mass flow of pebbles across the orbit of the considered core, 
regardless of how that flow had been reduced by the presence of additional cores further 
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out. We stress that the mass flow is related to the time at which a certain mass is achieved 
by a given body. Thus, reducing the mass flow over time (for instance because the cores 
outside of Jupiter intercept an increasing fraction of the overall flow) is just equivalent to 
stretching the x-axis of Fig. and in a non linear way. The masses of the inner solar 
system embryos as a function of the mass of the core of Jupiter would not change. 

However, one might wonder what would be the relative mass growth among the cores 
beyond the snowline. As we have seen, in the inner solar system the accretion rate of 
embryos grows with heliocentric distance. But the architecture of the Solar System, with 
Jupiter having accreted a more massive atmosphere than Saturn and Saturn having accreted 
more gas than Uranus and Neptune, suggests that the core of Jupiter grew the fastest, then 
Saturn’s, then those of the ice-giant planets. The growth of these cores has been already 
modeled in Lambrechts and Johansen (2014). They found indeed that the innermost core 
grows the fastest. The reason for this is twofold. The most important factor is that beyond 
the snowline the disk is flared, so that r] grows rapidly with distance (see Fig. [^). The 
second factor is that icy pebbles grow as they drift, so that their value of r increases with 
decreasing distance (see Lambrechts and Johansen, 2004, Fig. 2) which, in a moderately 
turbulent disk in which most of the accretion proceeds in a 3D mode, leads to more efficient 
accretion. Instead, within the snowline rj is almost constant and r decreases with decreasing 
r. Our calculations, assuming the //-function illustrated in Fig. and the r-function of 
Lambrechts and Johansen (2014) indeed confirm that, beyond 5 AU cores formed by pebble 
accretion have a growth rate decreasing with increasing heliocentric distance. 


4. Conclusions 

In this paper we have investigated which mode of planet growth can best explain the 


great dichotomy of the Solar System. This dichotomy stems from the realization that. 
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in order to form the giant planets, solid cores of approximately 10-20 Earth masses had 
to form in the outer solar system within the lifetime of the gaseous protoplanetary disk, 
whereas in the inner Solar System planetary embryos had a mass only of the order of that 
of Mars (as deduced from the short timescale of Mars accretion and from the mass of the 
embryo that was involved in the Moon-forming collision with our planet). How to get a 
ratio of ~ 100 between the masses of the protoplanets in the outer and the inner solar 
system respectively is not obvious. 

We have first considered the classic mode of protoplanet growth, that is ordered 
accretion of planetesimals aided by the so-called gravitational focusing factor. This leads 
to the process of runaway growth (Greenberg et ah, 1978), followed by that of oligarchic 
growth (Kokubo and Ida, 1998). We have shown in Sect. 2 that the runaway growth process 
is much faster in the inner disk than in the outer disk, in sharp contrast with the great 
dichotomy. Only if the process of oligarchic growth is brought to completion, with each 
embryo accreting basically all the solid material in its annulus of influence, one can expect 
more massive bodies further out. This is unlikely to happen, particularly in the giant planet 
region (Levison et ah, 2010, 2012). Nevertheless we find that, if oligarchic growth had 
reached completion everywhere, the mass of the embryos would have increased progressively 
with heliocentric distance. Thus, in order to have a body approaching IOMq beyond the 
snowline. Earth-mass embryos should have formed in the asteroid belt. Consequently, the 
resulting mass distribution in the Solar System would not have looked like anything that has 
been considered so far in successful models of giant planet or terrestrial planet formation. 

A second process for the formation of giant planet cores and planetary embryos is 
pebble-accretion (Lambrechts and Johansen, 2012). In that process, the most massive 
planetesimals accrete small objects (pebble-sized) as the latter drift through their orbits. 
Accretion of pebbles is very effective due to a combination of gravitational deflection 
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and gas-drag. We have shown that the process of pebble-accretion can explain the great 
dichotomy, provided that a nnmber of assnmptions hold trne. We have assnmed that icy 
pebbles are a few cm in size, according to the pebble-growth model of Lambrechts and 
Johansen (2014). When icy pebbles drift throngh the snowline the ice snblimates and 
we have assnmed that this causes a drop of a factor of 2 in the total mass flux and the 
release of a large number of silicate grains whose sizes are approximately 1 mm (i.e. about 
chondrule-size). With these assumptions we have been able to show that, starting from 
bodies of equal masses everywhere in the protoplanetary disk (specihcally, we have assumed 
a mass of 0.5 Lunar-masses), the embryos within the snowline grow to approximately 
Mars-mass in the time when the first core beyond the snowline grows to 2OM0. The 
major condition for this to happen is that there is enough turbulence in the disk so that 
the thickness of the layer of silicate grains is larger than the effective accretion radius for 
Mars-mass bodies. This happens if the parameter a characterizing the turbulence strength 
in the Shakura and Sunyaev (1973) prescription is > 10“"^. According to Stoll and Kley 
(2014) this condition is fulhlled due to the vertical shear disk instability. 

If this condition holds, the hnal mass distribution of the planetary embryos within 
the snowline moderately grows with heliocentric distance. Interestingly, this is the kind 
of distribution that reproduces best the final mass distribution of the terrestrial planets, 
including a small Mercury (Jacobson and Morbidelli, 2014). 

We kept the model simple on purpose, as a proof of concept. In particular, we have 
assumed that all the seeds of the future embryos and cores had initially the same masses; 
we have neglected the disk’s evolution; we have considered only one size for the pebbles in 
each region of the disk. Obviously, the results could easily change at a quantitative level 
by changing the assumptions. Thus, at this stage, in absence of more precise information 
on the original distribution of the largest planetesimals throughout the disk and the 
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sizes of pebbles and grains, we cannot be firmly predictive of what kind of protoplanets 
formed through the Solar System within the lifetime of the disk of gas. Nevertheless, we 
find interesting and intriguing that the pebble-accretion process, with a set of reasonable 
values for the parameters of the problem, allowed us to reproduce the great dichotomy 
of the Solar System, with Mars-mass embryos throughout the inner solar system and 
giant planet cores of IO- 2 OM 0 beyond the snowline. Given that the classical process 
of planetesimal-accretion did not allow us to come even close to the desired result, we 
conclude that the pebble-accretion process is more likely to explain correctly the formation 
of protoplanets throughout the Solar System. 

The processes of pebble growth and drift through a disk, as well as the pebble 
sublimation and splitting when passing through the snowline, should be universal. Thus we 
expect that any protoplanetary disk developed a sort of dichotomy between the masses of 
the protoplanets formed within and beyond the snowline, like the one we see (or, better, 
deduce) in our Solar System. If this is true, the great diversity observed among planetary 
systems would have to stem from the subsequent evolutions of these protoplanets. We 
believe that this is likely. For instance, if the cores beyond the snowline do not become 
giant planets, they would eventually migrate towards the star as the disk loses mass (Bitsch 
et ah, 2015). In this migration process the cores would scramble, disperse and deplete the 
system of planetary embryos, preventing the accretion of real terrestrial planets (Izidoro 
et ah, 2014); on the other hand the cores, moving into the inner system, would qualify 
for “super-Earths” or “hot-Neptunes” (probably low-density ones given their icy nature 
and their potential hability to accrete some primitive atmosphere), which are the most 
abundant planets observed around stars (Mayor et ah, 2011). If instead the cores become 
giant planets, their migration history becomes sensitive to their mutual interactions (Masset 
and Snellgrove, 2001). If only one core becomes giant planet then it should migrate inwards 
and this should strongly affect the evolution of the inner system of planetary embryos (Fogg 
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and Nelson, 2005; Raymond et al., 2006b). The key aspect of the Solar System evolntion 
was that Jnpiter avoided to migrate inwards of a few AU, probably dne to the presence 
of Satnrn (Masset and Snellgrove, 2001; Morbidelli and Crida, 2007; Walsh et ah, 2011). 
Thns, Jnpiter and Satnrn blocked the way to Uranns and Neptnne, preventing these cores 
to migrate into the inner Solar System as the disk cooled (Izidoro et al, 2015). This way, 
the giant planets remained in the onter Solar System, while the embryos snrvived in the 
inner Solar System, eventnally leading to the formation of a few terrestrial planets. This 
specihc evolntion preserved the vestige of the original dichotomy of protoplanetary masses, 
set in the primary pebble-accretion process. 
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